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一 类 2- 指 标 变 延迟 微分 代数 方程 BDF 方法 的 收敛 性 * 
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dS 要: 延迟 微分 代数 方程 经 常 出 现在 自动 控制 、 电 力 和 电路 分 析 、 多 体 动 力学 等 许多 实际 应 用 问题 中 ， 
目前 对 延迟 微分 代数 方程 数值 分 析 研究 主要 集中 于 线性 问题 和 1- 指 标 问 题 ; 对 高 指标 非 线性 延迟 
微分 代数 方程 数值 分 析 的 研究 较 困 难 ， 国 内 外 仅 有 少量 工作 且 大 多 为 常 延 迟 . 本 文 将 向 后 微分 公 
X (BDF) 应 用 于 求解 2- 指标 非 线 性 变 延迟 微分 代数 方程 ， 获 得 了 相应 的 收敛 性 结果 ， 并 通过 数 


值 试验 进行 了 验证 . 
关键 词 : 2- 指 标 微分 代数 方程 ， 向 后 微分 公式 ， 收 敛 性 ， 变 延迟 
分 类 号 : AMS(2000) 90C29 中 图 分 类 号 : 02214 文献 标识 码 : A 
1 引言 


目前 ， 对 延迟 微分 代数 方程 的 研究 思想 许多 借助 了 微分 代数 方程 的 理论 及 延迟 微分 方程 的 
理论 ， 因 此 常 微分 方程 ， 微 分 代数 方程 及 (中立 型 ) 延迟 微分 方程 中等 常用 的 理论 及 数值 方法 
被 许多 学 者 应 用 于 延迟 微分 代数 方程 及 其 数值 方法 . 但 是 ， 延 迟 微分 代数 方程 因 有 延迟 的 影 
响 ， 很 多 特性 是 微分 代数 方程 或 延迟 微分 方程 所 没有 的 ， 所 以 对 延迟 微分 代数 方程 的 研究 必须 
独立 进行 ， Petzoldlal 指出 “微分 代数 方程 不 是 常 微分 方程 ”， 同 样 ， 延 迟 微分 代数 方程 既 不 
是 微分 代数 方程 ， 也 不 是 延迟 微分 方程 . 

迄今 为 止 ， 国 内 外 对 延迟 微分 代数 方程 的 研究 集中 在 理论 分 析 及 数值 方法 的 稳定 性 和 误 
差分 析 上 . 在 理论 分 析 方 面 ， 重 点 研究 了 (中 立 型 ) 延迟 (积分 ) 微分 代数 方程 的 结构 分 析 及 
JERES. 在 数值 稳定 性 方面 ， 重 点 讨论 了 Rosenbrock、( 多 步 ) Runge-Kutta(RK) 方法 
及 线性 多 步 法 ( 含 6- 方 法 、BDF 方法 ) 关于 线性 (中 立 型 ) 常 延 迟 (积分 ) 微分 代数 方程 的 渐 近 稳 
定性 的 5,8.9 及 0- 方 法 ( 含 中 点 法 ) 关于 1- 指标 刚性 延迟 微分 代数 方程 的 稳定 性 [10, 蔬 ， 在 数值 误 
差分 析 方 面 ， 文 献 [12-14] 分 别 获得 了 RK 方法 、 线 性 多 步 方法 、 单 支 方法 及 块 方法 关于 1- 指 
标 、2- 指 标 常 延 迟 微分 代数 方程 的 误差 分 析 结 果 ， 文 献 [15] 获得 了 配置 方法 关于 1- 指 标 、2- 指 
标 变 延迟 半 显 式 非 线 性 微分 代数 系统 的 收敛 性 结果 ， 文 献 [16,17] 获得 了 单 支 方法 和 RK 方法 关 
于 1- 指 标 变 延 迟 微分 代数 方程 的 收敛 性 结果 . 

本 文 应 用 向 后 微分 公式 (BDF) 求解 一 类 2- 指 标 非 线 性 变 延 迟 微分 代数 方程 ， 获 得 了 该 方法 
的 收敛 性 结果 ， 并 通过 数值 试验 进行 了 验证 . 
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2 一 类 2- 指 标 变 延迟 微分 代数 方程 的 收敛 性 分 析 
考虑 一 类 2- 指标 变 延 迟 微分 代数 方程 (2-DDAEs) 
y (z) = f(y(z),y(x —7(z).z(z)). «€ [0,T), 
0 = g(y(z)), z € [0,T], (1) 
z(0)= zo, y(x) = p(z), z € [-7,0], 


这 里 延迟 函数 r(z) 可 导 ， 且 满足 0 < r(z) <《T 0 < r (£)<1, RA% f: R™ x R™ xR” 一 

有 BR g: R 一 Rn" 充分 光滑 且 各 阶 导 数 有 界 ， 初 值 函数 po(z) 连 续 ， 且 gv(y) fi(y, y (m 一 

T(z)),2) 可 逆 . 恒 设 问题 (1) 有 唯一 光滑 真 解 y(z), z(z). 在 下 文中 符号 外 . | 表示 欧 氏 范 数 ， 
把 带 p 阶 插值 公式 的 BDF 方法 应 用 于 求解 问题 (1)， 可 得 





k 
DD Qiynti = hf (Yntk, ył, Zn+k)» 0= g(Yn+k), (2) 
i—0 


XE Tnk 二 Tn 十 kh, n20. H 


p(Tntk — T(Tn+k)), Zak — T(Zn4k) < 0, 


q 
P» Qj (OL Wa kk mad Tn+k 一 T(Tn+k) > 0, 


j--u 
其 中 T(Zn4x) = (ma — ôn)h, u,q, ma € Z*, ôn € (0, 1), q+u =p, 4+1 < mn, Qj (ôn) 为 拉 格 
朗 日 插值 基 函 数 . 

考虑 (2) 的 扰动 问题 


k 
y» Aini = hf (ntk, È, ntk) + hó, 0= g(fuax) +9. (4) 


i=0 
定理 1 B yntk, zntk E (2) 确定， 而 扰动 加 +k, ntk EH (4) 确定， 起 始 值 满足 


Ynti — YTnti) = O(h), znti — 2(2n44) = O(h), 


gynti) =O), j=0,l ,kl Znti= En + jh, (5) 
及 起 始 值 扰动 满足 
ti — Vn. = O(h?), Bnti— 245 — O(h), 8=0(h), 8-—0(n), (6) 
则 对 h < ho, A 
lôntr — nli < C(IIYa — Yl + Rll + ol ), (7) 


k-1 
lisa = ensell S Z S lose Ges — | MES — Yali + Alal + el). (9) 
j=0 


这 里 
a A) (9) 
IY, = Yall = max (mg) (10) 
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证 明 $ 
k—1 ai . k—1 oi. 
7 — Br n= S n (11) 
调整 hn 和 6 比例， 使 得 (2) 和 (4) 分 别 变 为 
Ynt+k =N + hf(Yn+k, yË, Zn4k), 0- g(Yn+k), (12) 
Ûn+k = À + hf (ntk R, ntk) +h, 0 = gnt) +0. (13) 
(12) 中 第 二 式 可 改写 为 
1 
0 = g(yn+x) — 9(7) + 9(7) = f gy (n +T(yn+k — m)dr|ys-ex — n] + 9(7), (14) 
把 (12) 中 第 一 式 代 入 (14) 可 得 
1 
f gy (n+ (yn — m)dr f (Ynte Y, Zntk) + an) = 0, (15) 
同 理 ，(13) 中 第 二 式 改写 为 
1 
J sr arl Onn ++ + i6 ao) 
利用 函数 f, g 的 光滑 性 及 gy 亡 的 可 逆 性 ， 由 (15)-(16) 可 得 
| 和 一生 Ci (ld — ell + I — mL 19 — vl 
1 1 
十 | + FAN + ZNA) — s01), (17) 
& (12) 中 第 一 式 减 去 (13) 中 第 一 式 可 得 
ld — vnl € IA- nl t h|| f nos IR Ink) — F Untk Veo 2x) || 
< J- nli - AL (Ilse — Yarell + MIR — vell linn — nell), (18) 
这 里 工 为 右 端 函数 了 的 经 典 Lipschitz 常数 ， 把 (17) 代入 (18)， 化 简 得 
1 
linen — medi < Co (Il — nll + lk vele MISERE N01), hc C9 
把 (19) 代入 (17)， 化 简 得 
上 ts — zl < Co( II — nli + 88 — v + Iall 
1 1 "On 1 
HION 十 glAN- MD), hs zr (20) 
x 
k-1 y, : 
l= oll = || D es e| s eus - Yall; (1) 
了 一 0 
IMAA- M|] € lle) - m] + lâ — n0 C) 
k-—1 
< Cs( > lo) nss — e| IIS = YaloQ). — Q2) 
了 =-0 


338 I E 数 学 学 报 第 28 卷 


把 (21), (22) 代入 (19), (20)， 经 化 简 得 结论 (7), (8). 
推论 1 设 2-DDAEs(1) 满 足 gyfz 可 逆 ， 且 带 p 阶 插值 公式 (3) 的 BDF 方 法 (2) 是 p 阶 相 容 
的 ， 则 它 的 局 部 截断 误差 满足 


yk — y(zx) = O(hP*), zk — z(£ę) = O(hP). (23) 


A 令 定 理 1 中 n= 0 9; = ylz), 2; = z2(zj), j 2 0,1, ,k， 这 些 值 满 足 定 理 1 的 条 
件 ， 且 5 = O(h»), 6—0, WAR (3) 4 gè — yk = O(h?+1)， 应 用 定理 1 即 得 结论 . 
定理 2 设 2-DDAEs(1) 满 足 gyf 可 逆 ，k «6, 方法 (2) 起 始 值 满足 


yj — y(z) = O(h**1), j-0,1:-- ,k— 1, (24) 
则 带 p 阶 插值 公式 (3) 的 k 步 BDF 方 法 (2) 是 p = 大 阶 收敛 的 ， 即 
Vn — Y(n) = O(h?), zn —z(z,)—- O(hP)), zn =nh, n>k. (25) 


证 明 首先 研究 y 分 量 的 局 部 误差 积累 . 数值 解 序列 记 为 {y8, z0}, 考虑 多 步 解 序 列 (ll, 253, 
l = 1,2,.….， i& 16 (8 7h vi = y(zj z} = z(zj), j = 二 1 一 1,… ,1 十 kk 一 2. 首要 任务 是 估 
计 匈 一 1， 为 了 方便 起 见 ， 我 们 去 掉 上 标 ， 两 个 相 邻 多 步 解 分 别 记 为 (On, 2n} Win, Zn} 
为 了 能 利用 定理 1 的 结论 ， 固 定 3 个 充分 大 的 常数 Co, C1, Cs ( 文 后 注释 说 明了 这 3 个 常数 的 合 
理性 )， 使 


le; — 9m) € Coh, — les — z(zntj) € Gs 

ls; — tl < OSA, 0,1, k- 1. (26) 
引入 记号 AZn+k = Znk — Ên+k; AYn+j = Ün4-j — Dntj j 一 0, 1, E ,K, Ay? = gh - ih, AY, et 
(Ayfer AYZ, (Ayh), EZI PRS = 0, 9 = 0， 应 用 定理 1 可 得 


C k—1 : 
lAyntall S CCIIAYsII), MAznsell < $Ò 3 llo oe) Ael + MAYA). (27) 
j=0 


只 要 六 充分 小 时 ， 常 数 C 不 依赖 于 Co, Cu, C2 的 选取 ， 由 插值 公式 (3) 及 假设 (26)， 有 
lAyn+rll = OQ), ||Azn+rl] = O(h). (28) 


线性 化 多 步 公式 得 


k 
Moi = h[f (Bk do Enek) — f Unies 15 2n] 
i-o 


= hf (nr, 95,2 o5) A 2544  O(h|AY,I|), (29) 
id Qa; 一 (gy f2) 9y) Ôn, DÈ, n+), Py —-1-—Qma j = 0, l,- ks 易 证 


Q2; = Quo, P = PH Qn+jPntj = PuajQnej =0, Qn+j+1 = Qn; O(h), (31) 
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H Pa+k 左 乘 (29)， fz(gyfz)-! 左 乘 (30)， 可 得 


k 
5 Oi Ps ci yas T o( h||AY |l ), Qn+kAYn+k T O( || AY; | ). (32) 
i=0 


由 插值 公式 (3)， 延 迟 函 数 0 < r'(z) < 1 可 导 ， 假设 (26) 及 (27)， 有 
1 
5Avla = Auh  O(M]AY, IJ), (33) 
记 

T T 1 AATAT 
= ((BoekcaAtne-i) s (PaA) zh). ) ， 
T 
Va = ((Quis-iAge- 1) s (QnA), z (AW)” ) ， 


A AY, = Un + Vn (32) 式 可 分 别 化 为 


Usi = (4 8 DU, + O(h|Usl| + Alival), (34) 
Vai = (N 8 DV, + O(h|U | + Alival), (35) 
这 里 
-ap -oa -ao 0 0 .000 
1 0 0 0 1 000 
A- : |， N=|: ; 
0 1 0 0 100 
0 0 0 1 0 0 0 1 


其 中 oj = o;/ox. 由 于 p(6) = PE 类 似 于 文献 18j]， 可 选取 一 范 数 | . |， 使 
filAeI|«1 KA-T Ex igni |- b &&INer|zi AKE 


mo p O(h) qe | a5 

[Vasall Olh) 1+0(h) J V vl 

4t. (36) 中 系数 矩阵 成 对 角 阵 ， 易 得 
Er mL llUoll T 
IV. 0 X II Voll 


这 里 Ai = 1--O(h), A2 = 1--O(h), ZRET ELTHIESEAEISI EPI, WE 


.( 1 om 
A00 d f’ 
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化 简 得 
IU, ll < Č(|lUoll + hllvoll), 
[vall < C(hIUol| + lvoll), 
IY < Ce(IUsI + IVall) < Cr( Voll + llVoll )- 
HE Uo, Vo 由 y 变 量 的 起 始 值 的 误差 及 局 部 截断 误差 组 成 ， 满 足 
[Uol] < Hoh"*?, Voll < Hih?t!,  ||Ayol < Hoh***, 
由 (23), (24), (38)-(40) 及 (41)， 可 知 


1 Ayn| < Cahtt, lu (0o) Am; || < Coh?t!, 


n—k41 
llus — yzo) YO lyh — L1 < Cioh. 
{=0 
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(42) 


(43) 


由 于 zn 仅 依赖 于 内 sys t 3 Yn—-1>? 在 定理 1 "PU; zn y(zi), ĉi = z(zi), 0 = O(h?), 8 = 


0， 应 用 定理 1 的 结论 可 得 
k 
ln — z(2«)ll € — Soy (yzn)) yn; — Y(£n-5))|| + OC). 
320 


利用 (42) 和 (43)， 可 得 


n—k+1 


loy (u(25))9«-5 — ylen- = X Ilgai + ORP- — 9 3| 
1-0 
n—k-F1 


I^ 


因此 
zn. F^ z(zn)|| < Ch. 


(44) 


2- (llo. — yit] + o(h?P+1)) — O(hPt!), 


(45) 


注 1 一 般 来 说 ， 常 数 Cio, Cu 依赖 于 假设 (26) 中 的 Go, Cl， 我 们 严格 限制 步 长 h， 使 


得 Cigh?-! < Co，Ci1h?-1 < CG， 这 样 ， 数 值 解 就 不 会 与 假设 (26) 相 冲突 . 


3 ”数值 试验 
考虑 如 下 2- 指 标 变 延迟 微分 代数 初 值 问题 
yi (x) = 29312 + nE, 0<Z<2， 
ylz) = -£ — yis () — VZ, 0<z<2, 
0= yiya +4, 0<z<2， 


yı(0)=1, y2(0)=-—1, z(00 — 1, 
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问题 (46) 有 唯一 真 解 yi (m) = e77, yo(z) = -er-z，z(z) =e. RIK h, yi, yo, z 在 z = 2 处 
的 整体 误差 分 别 记 为 yerrl(h)，yerr2(h)，zerr(h). 通过 计算 可 以 得 到 相应 的 误差 阶 ， 分 别 记 


A pyl(h), py2(h), pz(h)， 其 中 
— In2. 


yerrl(h) 
分 别 用 带 线 性 插值 公式 的 二 步 二 阶 BDF 方 法 (BDF2) 和 带 二 次 插值 公式 的 三 步 三 阶 BDF 方 


_zerr(h) 
zerr(0.5h) 


yerr2(Rh) 


pyi(h) = In / 2, py2(h) 一 In eris / n? pz(h) — In 


yerr1(0.5A) 


法 (BDF3) 求解 初 值 问题 (46). 


数值 结果 如 表 1 和 表 2 所 示 ， 验 证 了 理论 结果 的 正确 性 . 


表 1: BDF2 的 数值 结果 
h yerri yerr2 zerr pyl py2 pz CPU 
0.1 0.3590E —2 0.4899E —2 0.7199E—2 1.9974 2.0283 2.0865 0.1590 
0.05 | 0.8991E —3  0.1201E —2  0.1695E —2 2.0031 1.9750 1.9924 0.3167 
0.025 | 0.2243E —3  Á0.30558E —3  0.4260E —3 2.0019 1.9976 2.0122 0.6085 
0.0125  0.0560E —-3  0.07658E —3  0.1056E — 3 1.1973 
表 2: BDF3 的 数值 结果 
h yerrl yerr2 zerr pyl py2 pz CPU 
0.1 0.5100 -3 0.700E -3 40.1100E —2 2.9521 2.9436 3.0179 0.2179 
0.05 | 0.6590E —4 0.91008 —4  0.1358E —3 2.9684 3.0019 3.0541 0.4064 
0.005 0.8420E —5  0.1136E —4  0.1685E —4 3.0076 2.9501 3.0589 0.7932 
0.0025 0.1047E —5  0.1470E —5  0.1962E — 5 1.4571 
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Abstract: Delay differential-algebraic equations often arise in automatic control, power and circuit 
analysis, multi-body dynamics, etc. The current researches on numerical analysis for delay-differential- 
algebraic equations are mainly focused on linear problems and l-index problems. It is difficult to do 
numerical analysis for high-index nonlinear delay-differential-algebraic equations, and there are only a 
few results for this kind of problems, furthermore, most of them are about constant-delay problems. The 
backward differentiation formulas are applied in this paper to index-2 nonlinear differential-algebraic 
equations with variable delay. The corresponding convergence results are obtained and confirmed by 
some numerical examples. 
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